ro = 90;
col = 120;
% new C, lb, ub, b
c5=c2(1:col,1);
lbounds5=lbounds2(1:col,1);
ubounds5=ubounds2(1:col,1);
AT5 = AT2(1:ro,1:col);
bT5 = bT(1:ro,1);

lbConstraintA4=lbConstraintA(1:col,1:col);
lbConstraintb4=lbConstraintb(1:col,1);

tic;
xlswrite('D6CUBE-90x120.xlsx', c5, 'c');
xlswrite('D6CUBE-90x120.xlsx', AT5, 'A(<=)');
xlswrite('D6CUBE-90x120.xlsx', bT5, 'b(<=)');
xlswrite('D6CUBE-90x120.xlsx', lbConstraintA4, 'A(>=)');
xlswrite('D6CUBE-90x120.xlsx', lbConstraintb4, 'b(>=)');
xlswrite('D6CUBE-90x120.xlsx', Aeq, 'Aeq');
xlswrite('D6CUBE-90x120.xlsx', beq, 'beq');
toc

clc;
tic
% [x3,fval,exitflag] = linprog(c,A,b,[],[],lbounds,ubounds);
% [x3,fval,exitflag] = linprog(c2,AT3,b3,[],[],lbounds2,ubounds3)
% [x3,fval,exitflag] = linprog(c5,AT5,bT5,[],[],lbounds5,ubounds5)
[x3,fval,exitflag] = linprog(c5,AT5,bT5,[],[],lbounds5)
toc